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ABSTRACT 

We present the results of a MIPS-24^m study of the Brightest Cluster Galaxies (BCGs) of 535 high- 
redshift galaxy clusters. The clusters are drawn from the Spitzer Adaptation of the Red-Sequence 
Cluster Survey (SpARCS), which effectively provides a sample selected on total stellar mass, over 
0.2 < z < 1.8 within the S'piteer Wide-Area Infrared Extragalactic (SWIRE) Survey fields. 20%, or 
106 clusters have spectroscopically confirmed redshifts, and the rest have redshifts estimated from the 
color of their red sequence. A comparison with the public SWIRE images detects 125 individual BCGs 
at 24^m > 100/rJy, or 23%. The luminosity-limited detection rate of BCGs in similar richness clusters 
(A"gai > 12) increases rapidly with redshift. Above z ~ 1, an average of ^ 20% of the sample have 
24/rm-inferred infrared luminosities of L/r > 10^^ Lq, while the fraction below z ~ 1 exhibiting such 
luminosities is < 1 %. The 5pzteer-lRAC colors indicate the bulk of the 24/rm-detected population 
is predominantly powered by star formation, with only 7/125 galaxies lying within the color region 
inhabited by Active Galactic Nuclei (AGN). Simple arguments limit the star-formation activity to 
several hundred million years and this may therefore be indicative of the timescale for AGN feedback 
to halt the star formation. Below redshift z ~ 1 there is not enough star formation to significantly 
contribute to the overall stellar mass of the BCG population, and therefore BCG growth is likely 
dominated by dry-mergers. Above z ^ 1, however, the inferred star formation would double the 
stellar mass of the BCGs and is comparable to the mass assembly predicted by simulations through 
dry mergers. We cannot yet constrain the process driving the star formation for the overall sample, 
though a single object studied in detail is consistent with a gas-rich merger. 

Subject headings: galaxies: clusters: general, galaxies: evolution, galaxies: high-redshift, galaxies: 
starburst 


1. INTRODUCTION 

At the centre of most local galaxy clusters lies a 
single massive galaxy. These so-called Brightest Clus¬ 
ter Galaxies (BCGs) are the most massive galaxies in 
the universe today. They do not appear, however, to 
simply be the extremes of the local galaxy mass func¬ 
tion, but are separate beasts, with luminosities, metal- 
licities and surface brightness profiles t hat are unique 
and consistent across the populatio n (lOemlerl 1 1 9761 

Tre i naine fc Hichstonel 119771: iDressled Il978t ~hombertl 

1986h . It seems likely that the formation of BCGs is tied 
to the overall growth of their parent clusters, through the 
physics of gas cooling in the most massive dark matter 
halos and their hierarchical accretion of the general field 
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galaxy population. 

Nevertheless, we know rather little about the for- 
mation histor i es of BCGs. Recent measurements by 
iLidman et al.l (I2012D show that BCGs h ave increased 
their stellar mass by 2 x since z ^ 1 Ic.f iCollins et al.l 
120091 : l^ott et al.l 1201 Ih . These authors find the growth 
is driven by dry accretion of sa tellite galaxies, in 
line with Semi-A nalytic Models (I Bower et al.l 120061 : 
iDe Lucia fc Blaizotl 120071 : lGuoll2010ll . The accreted sys¬ 
tems are gas-poor and no new star formation is in¬ 
duced during the merger process. In this picture most 
of the stellar mass is formed at very high redshifts 
(z ~ 5) within multiple galaxies through low levels of 
star formation, and the BCG slowly acquires its iden¬ 
tity through the conglomeration of previously assembled 
pieces. On the other hand, several moderate redshift 
BCGs (z < 0.6) exhibit signs of star formation (lO-IOOO 
M|^yr~^) and/or contain large amounts of molecu l ar gas 
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2007; 

O’Dea 

201C 

Donahue et al.l 120111 IRawle et al.l 

2012; 

McDonald et al.ll20I2r). indicating the scenario is not so 


simple. In many cases the star formation rates are corre¬ 
lated with the gas cooling time, and indeed star forma¬ 
tion has been primarily seen in cool-core clusters. It is 
not yet clear how important this star formation is to the 
overall growth of BCGs as a population, nor how it might 
relate to similar processes at very high redshift that are 
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posited by the models. 

Star formation fed by rapid gas deposition and contin¬ 
ual dry galaxy accretion are two very different scenarios 
of BCG evolution, and each have implications for the 
physics of halo collapse and growth. They are not mutu¬ 
ally exclusive since both are transient events regulated by 
different processes and may occur at different times in the 
history of a galaxy. BCG galaxy accretion is an ongoing, 
though sporadic, phenomenon, fed by the continual in¬ 
fall of field galaxies into the cluster potential, which are 
stripped of much of their interstellar medium (ISM) dur¬ 
ing their descent into the centre of the cluster. Cooling 
flows, and the resulting central starburst, likely have de- 
hned duty-cycles that are go verned by still unclear heat¬ 
ing and cooling mechanisms (iMcNarnara fc Nulsenll2007l : 
iRaffertv et al.ll200^ IVoit et al.ll2008[l 

Studies of galaxy clusters and their central galaxies 
suffer from observational bias and, until recently, small 
sample sizes. Much of the work on central star forma¬ 
tion has been carried out on X-ray selected clusters which 
may preferentially select cooling flow clusters. Quite re¬ 
cently, large samples of Sunyaev-Zeldovich clusters have 
become availab le through facilities such as the South Pole 
Telescope fe.g. [McDonald et al.ll2013^ . that also select 
clusters through observations (though indirect) of their 
intercluster medium. Here we investigate the observed 
IR properties of BCGs using a large sample of clusters 
from the Spitzer Adaptation of the Red-Sequence Clus¬ 
ter Survey (SpARCS), an optical/NIR galaxy-selected 
cluster survey , designed to discover clu sters to z ^ 2 
(|Muzzin et ai.l l2009t IWilson et ^ l2009tl . While an im¬ 
portant expansion into new parameter space, this ap¬ 
proach also introduces new limitations and challenges. 
With a large sample of galaxy clusters considerably less 
is known about each individual system and the conclu¬ 
sions are generally statistical. It is necessary to automate 
some aspects of the analysis which increases the risk of 
sample contamination. Thus, we present the results of 
this work with the reminder that ongoing studies, such 
as high-resolution imaging and additional spectroscopy, 
are required to better calibrate some of the assumptions 
made here, and relate the findings to those of X-ray and 
SZ-selected samples. 

The paper is outlined as follows. In §2 we introduce 
the SpARCS cluster sample. In §3 we outline our BCG 
identification algorithm and the IR analysis using the 
SWIRE data. In §4 we discuss the origin of the IR emis¬ 
sion and the change in the IR properties of BCGs with 
redshift. §5 finishes with a discussion of the implications 
of §4 and we present our conclusions in §6. Standard cos¬ 
mology (Ho = 70 km/s/Mpc; ^matter = 0.3; Ha = 0.7) 
is adopted throughout. 

2. DATA AND OBSERVATIONS 

2.1. The Spitzer Adaptation of the Red-Sequence 
Cluster Survey 

The Spitzer Adaptation of the Red-sequence Cluster 
Survey (SpARCS) is a deep z'-band imaging survey cov¬ 
ering the Spitzer SWIRE Legacy fields, designed to as¬ 
semble a large homog eneously-selected sample of mas- 


sive clusters to z ~2 (IMuzzin et al.l 12000 

; IWilson et al.l 

120091 IDemarco et al.l 120101 IMuzzin et al.l 

201211. It em- 


ploys an infrared adaptation of the cluster red-sequence 


method (|Gladders fc Yeel [2005D : using the z' — 3.6/xm 
color, which spans the 4000A-break at z > 1, it locates 
over-densities of red-sequence galaxies that trace massive 
galaxy clusters. This simultaneously provides a reliable 
redshift estimate for the cluster, through the color fit to 
the red-sequence, discussed in more detail in the Section 
3. 

SpARCS uses the public SWIRE images and catalogs 
which provide photometry measurements at 3.6, 4.5, 5.8 
and 8.0 /rm from the IRAC instrument and 24/rm with 
the MIPS camera (Londsdale et al. 2004; Shupe et ah, 
private communication). SWIRE contains northern and 
southern sky coverage; however to ensure uniform an¬ 
cillary data (described below) we limit this work to the 
34 square degrees of the northern fields (ELAIS-Nl/2, 
XMM-LSS, Lockman), and further require uniform cov¬ 
erage in IRAC 3.6/im and 4.5/im. 

The parent sample is limited by the detection reliabil¬ 
ity of the red-sequence over-density and the measured 
NIR richness. Galaxy clusters must have a “flux” of 
at least 4 in over-densit y in the detection map - a cri¬ 
teria similar to that of iHildebrant et'M] (1201 Ih - and 
have a measured richness of Ngai >12. Ngai is de¬ 
fined as the number of background subtracted galaxies 
within an aperture of 500kpc above a limit of 
-|-1. These criteria ensure a richness-limited sample of 
high signal-to-noise detections (Figure [IJ. To place the 
richness in physical context we n ote that Ng gi = 12 cor- 
responds to M^nn ^ lxl0^^Mf7^ (iWen et alll2ni(Tl 120121 : 
iCapozzi et al.l l2012t lAndreon fc CongdonI 1201 J T The 
sample of 559 galaxies is further culled as discussed in 
Section 3. 

3. ANALYSIS 
3.1. BCG Identification 

Given the large number of clusters (a final sample of 
535, 106 confirmed and 429 candidates) within our sam¬ 
ple, we automate the BCG identification algorithm. We 
select the brightest galaxy in the IRAC 3.6/xm chan¬ 
nel whose color is within ±0.5mag of the predicted red- 
sequence. The location red- sequence is defined by the 
z' — 3.6R m color mode l s of Bruzual &: Charlotl (120031 1 
following iMuzzin et al.l (12009 ). Using clusters in com¬ 
mon we check this against the s maller, but more det ailed, 
study of high-redshift BCGs of (jLidman et al.ll20O) . and 
find that they select the same objects. This method 
is, however, susceptible to contamination by foreground 
dusty galaxies which are likely to be bright at 24/im and 
may have colors sufficiently reddened so as to be con¬ 
sistent with the red-sequence at higher redshift. Since 
the goal of this work is to study the infrared properties 
of the BCGs it is imperative that the frequency of such 
interlopers be quantified, and the misidentified BCGs re¬ 
moved from the analysis. 

We identify candidate interlopers through visual in¬ 
spection of the galaxy morphology. We use the Sloan 
Digital Sky Survey (SDSS) as it provides uniform imag¬ 
ing over all four fields. We identify galaxies with obvious 
spiral structure and remove these from the analysis. We 
perform a secondary check using the the deeper imaging 
of the Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS) in the XMM-LSS field, but do not identify 
any additional interlopers. In total, 24/559 galaxies are 
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redshift 



Fig. 1. — Top: The redshift distribution of the SpARCS cluster 
sample used here. Bottomt: The richness distribution (parameter¬ 
ized by Ngai) of the z < 1 clusters (black) and the z > 1 clusters 
(orange). Although fewer clusters are found at higher redshift, 
there is no significant difference in the richness distribution be¬ 
tween the two redshift bins. 


removed leaving a final sample of 535 BCGs. 

We then compare the spectroscopic redshift of the 
candidate BCG when available with the red-sequence- 
estimated redshift of the cluster. This comparison is 
shown in Figure H] and includes spectrosc opic 106 red- 
shifts assembled from SPSS , our own work (iMuzzin et abl 
l2012tl. the OzDES survey (iFange et al.ll2015D . PRIMUS 
( Coii_et_alJj20n|_lCool_et_^ 2013) , and the literatur e 


( Rowan-Robinson et al.ll20'08 , 12013 : (Wenger et al.l[2000l) . 

A systematic offset of z(RS) - 0.15 has been applied to 
the entire sample. This offset was empirically determined 
from these 106 redshifts to correct for differences between 
the model and measured colours of the red sequence. We 
are, as of yet, unsure of the source of this offset, but this 
will be further explored in Muzzin, A. et ah, in prep. As 
can be seen in Figure [5] there is excellent agreement be¬ 
tween redshifts and no additional galaxies are removed 
based on redshift disagreements. Although there are sev¬ 
eral BCGs whose RS redshift estimates are low compared 
to the overall scatter within the population, there is only 
one catastrophic outlier (|Az| > 0.5). This small offset 
may reflects a tendency of the RS method to occasionally 



RS-z 


Fig. 2.— A comparison of the cluster redshift, estimated through 
the location of the red-sequence, and the BCG spectroscopic 
redshift (open circles) or cluster spectroscopic (open stars) red- 


20081: ICoil et alj 120111: IMuzzin et alj 

20121: IGool et al.l 120131: 

Rowan-Hobinson et al.ll2013l: Ik'anee et al. 

|2015|). The solid line de- 


notes a 1:1 correspondence. The BCG catalog has been cleaned 
by-eye to remove all obvious low-redshift galaxy mis-identifications. 
Solid orange points show BCGs with 24/im counterparts (Section 
3.3), and these do not have a systematically different redshift offset 
than the 24/rm non-detected population. The overall scatter in red 
sequence redshift is cr = 0.10, using 106 spectroscopic redshifts. 


overestimate the redshift, rather than a misidentification 
of the BCG. After all interlopers are removed, we are left 
with a final sample of 535 BGGs. The richness of their 
parent clusters span Ngai =12 — 40 at all redshifts and 
the redshift distribution is shown in Figure [T] 


3.2. Caveats on Identification of the Brightest Cluster 

Calaxy 

This method of BGG identification has two additional 
biases. Firstly, it is biased against BGGs with unusual 
colors - i.e. colors that are inconsistent with the red- 
sequence. Thus, if a BCG were exceptionally blue or 
red, for example due to star-formation (unobscured or 
obscured), it would be passed over by the algorithm, and 
instead the second brightest red-sequence galaxy in the 
cluster would be chosen (maintaining the redshift agree¬ 
ment but technically selecting the wrong galaxy). In con¬ 
trast to foreground contamination, which adds spurious 
star forming galaxies to the sample, this effect may re¬ 
move real star-forming BCGs from the sample. 

Secondly, the sample could be contaminated by late- 
stage mergers that are unresolved in the IRAG-3.6^m 
imaging (with a Point Response Function of 1.8"). Such 
systems, which are actually two galaxies, but still cluster 
members, may appear as the brightest galaxy within the 
cluster when their light is combined. This latter effect 
would be more important at higher redshift where the 
apparent separation is smallest. We see no signs of this 
in the optical imaging, but higher resolution observations 
of a significant fraction of the sample are required to 
confirm this. 

Aside from the practical difficulties in identifying the 
BGG of a distant galaxy cluster, we must also consider 
the evo lving definition of the BC G itself. N-Body and 
SAMS (|De Lucia fc: Blaiz'otll2007D indicate that BGGs do 
not exist as a single identifiable object above z ~ 0.7, and 
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TABLE 1 

Information on cluster sample and observations 


SWIRE Eield 

Number 
of BCGs^ 
in sample 

Number 
of BCGs 

with S 24 ^.m > 100/iJy 

expected 

spurious 

matches 

XMM-LSS 

124 

29 

2 

Lockman 

202 

43 

2 

ELAIS-Nl 

145 

36 

2 

ELAIS-N2 

64 

17 

1 

total 

535 

125 

7 


^ With the cluster richness > 12. 


this is supported b y some, but not all, observations of dis - 
tant clusters (e.g. ISantos et al.l 120091 : IWebb et al.ll2015n . 
In these models, where the BCG assembles through dry 
mergers and galaxy cannibalism, the stellar mass which 
will eventually form the BCG is distributed among more 
than one galaxy. Even so, these galaxies will still be 
among the most massive systems in the cluster and our 
algorithm would select the brightest of them. Thus, our 
approach may be more accurately described as identify¬ 
ing one of the progenitors of the BCG galaxy. An investi¬ 
gation of the uniqueness of BCG galaxies at high redshift 
is possible with these data, but beyond the scope of this 
paper and will be the focus of later work by our group. 

3 . 3 . 24lJ>'m-detected BCGs 

The BCG catalog was cross-referenced with the 
SWIRE IRAG-24/im imaging catalogs, with a MIPS- 
BGG search radius of 3", to account for the uncertainty 
in the MIPS positions due to the 6" Point-Response- 
Function. Table I lists the number of 24/xm-detected 
BCGs in each field, to the depth of the SWIRE images 
(~ ISO/rJy), with a total of 125 detected systems, or 
^23% of the galaxies. Table 2 provides the coordinates, 
redshifts, and 24/im flux measurements of these systems. 

3 . 3 . 1 . Spurious Alignments 

We estimate the number of spurious spatial coinci¬ 
dences through simple Monte-Garlo simulations. Using 
the 24^m catalog we repeat our BGG matching, but 
with random positions in place of the BGG locations. 
The number of expected spurious detections are listed 
for each field in Table 1 and indicate a spurious BCG- 
24/Lm contamination rate of ~ 1% - much lower than our 
detection rate. We therefore do not attempt to correct 
for it in further analysis. 

3 . 3 . 2 . Physical Misidentifications: Gravitationally Sensed 
Systems and Close Neighbours 

We classify as physical misidentification cases where 
the IR emission is not associated with the BCG itself, 
but with an object which is physically related to the 
BGG. Here we consider gravitational lensing and close 
neighbours. 

Strong gravitational lensing from the galaxy cluster po¬ 
tential and/or the central massive galaxy can magnify 
background IR galaxies along the line of sight to the clus¬ 
ter centers. If the offset between the lensed object and 
the BCG is less than ^3" the IR flux will be confused 
within the Spitzer-MIPS beam and erroneously assigned 
to the BCG through the matching procedure. 



Fig. 3. — The HST-ACS-606W image one of the BCGs in this 
sample (from the CANDELS UDS public data) that is clearly lens¬ 
ing a background galaxy. In red we show the MIPS flux contours 
starting at 0.05 MJy/sr and increasing in 0.05 steps. Given the 
beam size of MIPS (6^^ FWHM shown) the source of the 24pm 
emission is ambiguous. 

To search for gravitational lenses within the SpARGS- 
BGG sample we turn to archival or public HST imaging. 
Very few (10) of the sample presented here have been 
covered in past surveys, but that number increases if we 
include the less rich ( but still meeting the significance 
criterion of Section 2.1) systems in SpARGS. In doing 
so we find one example of a tight (~2") lensing system 
in one cluster (Figure [3]). If this rate (~10%) is repre¬ 
sentative of the entire sample, lensing would indeed be 
an important source of contamination. However, this 
object is a previously known system that was first dis¬ 
covere d as a strong galaxy-galaxy lens by iGeach et al.l 
(|2007tl and later imaged in UDS field of the CAND ELS 
survey (|Grogin et al.l 1201 ll : iKoekemoer et al.l 1201111 . Its 
HST coverage may not be by chance, and therefore this 
rate may not be representative of the lensing rate. 

To gauge if this lensing rate is at all reasonable, given 
what we know about the system s, we ran lensing sim - 
ulations following the method of iHezaveh et al.l (j2012fl . 
We find the small radius lensing rates of only ^ 1%, 
which are too low to account for the IR detection rates 
of 23%, and moreover, the rate does not increase with 
redshift. While the lensing rates can be fine-tuned by 
altering the simulations (for example by tuning the el- 
lipticity or BCG-cluster alignments), we conclude that 
lensing is not a major contaminant of our sample. 

A second possible source of contamination that is phys¬ 
ically interesting is the presence of close neighbours to 
the BCG. By this we refer to physical neighbours, rather 
than line-of-sight associations which would be accounted 
for in our Monte Garlo simulations of spurious detections. 
These galaxies need to be close enough to be confused 
within a single MIPS 24/im PRF, but separated into two 
objects in IRAC 3.6/im. At such close separation the ob¬ 
jects will likely be interacting with the BGGs and thus 
since the two galaxies in a sense form a single system 
that includes the BCG, they modify but do not negate 
the interpretation of our results. 

To assess this issue, we have compared the frequency of 
close galaxy neighbours in the 24/im-detected and unde¬ 
tected samples using the full IRAC SWIRE galaxy cat- 
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alog and see no difference between the two, however a 
proper assessment of this effect requires higher resolution 
imaging so that very close neighbours can be revealed. 

4. RESULTS 

4.1. MIR Diagnostics of the Energy Source: AGN 
versus Star Formation 

In Figure |4] we show the 4-channel IRAC col¬ 
ors for the 24^m-detected and the 24^m- undetected 
BCGs. The IRAC colour plot is frequently invoked 
as a tool to discriminate between AGN (continuum- 
dominated), star-forming (PAH-domi nated) and passive 
(stellar-dominated) gala,xY systems (iSaiina et al.l l2005t 
Lacv et aH 12004 I2007t IStern et all I2005I iDonlev et ^ 
201^ . In Figure 4 we show the divisions suggested by 
Saiina et al.l (1200,^ (1: low-z PAH; 2 and 3: mid/high- 
z PAH and stellar continuum; 4: high-z PAH and all-z 
AGN) to statistically separate these populations as well 
as the more recent, an d tighter, AGN selection criteria 
of iDonlev et al.l (120121 1. 

The IR-faint BGGs lie almost entirely within regions 2 
and 3, indicating that they are truly passive systems, or 
higher redshift star forming galaxies with star formation 
rates too low to be detected in SWIRE. The low red- 
shift star-forming region is relatively empty, containing 
only one IR- undetected and seven IR-detected systems. 
The cluster redshifts for these systems are z < 0.6, in 
agreement with the expected low redshifts. The bulk 
(^75%) of the IR-detected objects lie within regions 1, 
2 and 3 with higher redshift systems preferring region 3; 
in agreement with the interpretation that their IR emis¬ 
sion is star-formation dominated. Region 4 contains only 
(bar one object) IR- detected systems. In this region 
we find a combination of high-redshift star-formers and 
AGN. Indeed, many of these systems form an extension 
of the colors of regions 2 and 3, suggesting they are also 
dominated by starbursts. 

Subsequent analysis identify as AGN only those 7 
systems which lie within the updated constraints of 
IDonlev et al.l (j2012( l. but otherwise include BCGs in re¬ 
gion 4. An upper-limit to the contamination from AGN 
to the MIPS detections can be obtained by considering 
all galaxies in region 4 as AGN. 

4.2. Mid-Infrared Luminosities of BCGs 

We convert the observed 24/im fluxes to inferred total 
IR lum inosities following the models of iGharv fc Elba^ 
(|200Ill . This assumes that all of the 24/rm flux is emit¬ 
ted by the BGG itself (that is, no gravitational leasing) 
and is entirely due to star formation with no AGN con¬ 
tamination. We therefore highlight, when relevant, those 
galaxies which may contain AGN, as identified in the pre¬ 
vious section (Donley et al. 2012) 

In Figure [S] we show the infrared luminosity of the 
24/im-detected BGGs with redshift and we list the val¬ 
ues in Table 2. We include fo r refer ence the Herschel 
detected BCGs of iRawle et al.l (j2012fl and the extreme 
star forming BGG in the Phoenix cluster reported by 
iMcDonald et al.l ()2012f ). The Phoenix BCG is an out¬ 
lier at z = 0.6, whereas the IR luminosities of the Rawle 
et al. sample are comparable to ours, with a few excep¬ 
tions. Note however, that the Rawle et ah clusters are 
X-ray selected and located at redshift z < 1.0, whereas 



-0.4 -0.2 0.0 0.2 0.4 0.6 


Fig. 4.— The IRAC color-color plot for the SpARCS BCGs. Blue 
points correspond to 24/im detected BCGs, and orange to 24/im 
non-detected BCGs. The increasing sizes of the points corresponds 
to increasing redshift. We also show 24/im detected BCGs that 
have limits on their 5.8/im or 8.0/im flux, but for clarity only in¬ 
clude the infrared-faint BC Gs that have detecti ons in all four IRAC 
channels. Overlaid are the ISaiina et all II2Q05I') SED template sec¬ 
tions (dotted divisions; 1: low- 2 : PAH; 2 and 3: mid/high- 2 : PAH 
and stellar continuum; 4 : high - 2 : PAH and all- 2 : AGN) and the 
more recent fPonlev et al.l 1)20121 ^ AGN wedge (solid line). 



redshift 

Fig. 5. — The 24/im-inferred inf rared luminosity (assuming 
the models of ICharv &: Elba3 <120011) .) Solid points refer to the 
SpARCS cluster sample (this work), with the orange colour denot¬ 
ing possible AGN, as identified by the Donley et al.(2012) region in 
Figure[4| The error bars include the 24/im flux uncertainties and a 
diO.l redshift scat ter. The open circl es correspond to the Herschel 
detected BCGs of IRawle et al.l 1120121) . although we use their 24/im 
estimated Lir for consistency. The open star denotes the extreme 
star-forming BCG reported by HMcDonald et aDl2012l) . discovered 
in the SPT survey. The solid star corresponds to the BCG with 
a strongly lensed arc (Section 3.3.2); here we assume the source 
of the 24/im flux is the BCG itself. Note that this object is not 
actually part of our study as it has Ng^i < 12, but is shown for ref¬ 
erence. The orange solid line shows the luminosity depth discussed 
in Section 4.2, that is tied to an evolving IR luminosity function, 
and the solid black line the approximate luminosity depth of the 
SWIRE 24/im survey. 

our sample is effectively stellar-mass selected to higher 
redshift, and therefore a direct comparison of the two 
samples is not straightforward. 

There is a consistent rise in Figure [5] in the inferred in¬ 
frared luminosities of BCGs with similar cluster richness 
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Fig. 6.— Top: The fraction of 24^m detected BCGs in different 
redshift bins. Blue points denote BCGs with inferred infrared lu¬ 
minosities of >Ljii 10^^ Lq and orange points denote the fraction 
of BCGs above the evolving luminosity cut shown by the orange 
line in Figure [5] Also shown (green) are the fractions of BCGs 
detected at any 24/rm flux above the detection limit of the SWIRE 
survey. Bottom: Same as above, but with possible AGN contami¬ 
nated galaxies removed from the analysis. For this plot we define 
as AGN those galaxies lying in region 4 of Figure 4. 

out to z ^ 2. Given that the infrared luminosity func¬ 
tion of field galaxies evolves toward higher luminosities 
and densities over this redshift range, this evolution may 
simply be a reflection of this general trend. In Figure [5] 
we show an evolving infrared luminosity limit following 
the Lir ~(1-|-z)^-^ relation of Le Floc’h et al. (2005) and 
scaled to our highest luminosity depth = 1 x 10^^ 

Lq at z = 1.8). This defines the luminosity limit as 
the same fraction of an evolving T/r, at every redshift, 
where, for reference ~ 1.6x10^° Lq at z ~ 1. 

There is still an increase in the absolute number of 
objects above this line, however this should not be inter¬ 
preted as a luminosity evolution; the absolute numbers 
of BCGs at a given luminosity will depend on the vol¬ 
ume probed (which increases by ^ 3 from 0 < z < 1 to 
1 < z < 2), as well as the cluster number density and 
selection function. 

In Figure [5] we show the fraction of BCGs that are 
detected above a given luminosity limit which is more 
likely to trace a real change in clusters with time. In 
the top panel of Figure |6] we show the fraction of all 


BCGs detected at different luminosity/flux limits. We 
include the fraction of BCGs detected at any flux above 
the depth of SWIRE (^lOO/rJy); above Lir > 10^^ Lq] 
and above the field related evolving infrared limit shown 
in Figure m In all cases the fraction of IR bright BCGs 
increases significantly beyond z ^I, for similar richness 
clusters. In particular, the field-corrected evolution in¬ 
creases from ~ 5% to ~ 30%, from the lowest to highest 
redshift bins. This plot includes all IR detected systems, 
even those with possible AGN contamination. In the 
lower panel we show the trend after removing all galax¬ 
ies within the Sajina et al. (2005) region 4, which as 
explained in Section 4, we take to be the highest level 
of AGN contamination. Using the Donley et al. region 
provides a much lower contamination rate. While the 
evolution is reduced by roughly a factor of two (for the 
Sajina et al. contamination rate), it remains significant, 
in part because the AGN are scattered throughout most 
of the redshift bins. 

We note, finally, a systematic effect that may 
be present in t his analysis. S everal studies have 
shown that the iGharv fc Elbad (j2001f ) methodology 
over-estimate s the infrared luminosity above a redshift 
of z ^ 1.5 (iMurnhy et ^ 120091 iNordon et al.l 120121 : 
iRodighiero et al.l I2010D bv roughly a factor of 5, due 
to large PAH equivalent widths at these redshifts. This 
redshift transition is the approximate location where we 
see the largest change in the fraction of 24/xm-detected 
BGGs. It may be that, due to this effect, we are sen¬ 
sitive to a lower luminosity limit beyond z ~1.5, which 
would in turn lead to a larger fraction of detected galax¬ 
ies. If we applied a downward correction to the L/r in 
Figure 5, this would indeed flatten the high-redshift tail 
and bring those estimations more in-line with the mea¬ 
surements below z ~1.5, but it would not change our re¬ 
sults qualitatively as these systems would still remain at 
L/r >1O^^L0. Moreover, the infrared luminosity of the 
most luminous object in our sample, SpARCSI049-|-56 
at z = 1.7, has been studied in detail in Webb et al. 
(2015), where we determine its luminosity using 6 in¬ 
frared measurements and two limits. We determine the 
L/r to be 6.6±0.9xl0^^ Lq, compared to an L/r to be 
1.0±0.3xl0^^ Lq determined from the 24/im measure¬ 
ment alone. These estimates are within Ict of each other 
indicating, at least for this particular object, there is no 
strong over-estimate of the flux. 

4.3. Star Formation Rates of the BCGs 

We can scale the infrared luminosities determined in 
the previo us section to star f ormat ion rates following the 
relation of iKennicutt et all (|1998D . Based on Figure |4] 
we assume no AGN contribution to the infrared flux, 
which would systematically reduce the estimates. We 
show the SFRs for each BGG with redshift in Figure [71 
Uncertainties include the photometric uncertainty in the 
MIPS 24/im flux and the redshift uncertainty of AziO.l. 
A wide range in star formation rates is seen, from ^1 at 
low redshift to ^1000 MQyi~^ for the highest redshift 
BGGs. 

Using stellar mass measurements for the BCGs we can 
further convert the star formation rates to star forma¬ 
tion rates per unit stellar mass, or specific star formation 
rates (sSFRs). Given the uniformity of the BCG optical 
colours we adopt a simple methodology to measure stel- 
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Fig. 7. — The inferred star formation rates of the BCGs with 
redshift. The orange points correspo nd to those wi t h IRA Q colours 
consistent with AGN following the IDonlev et al.l 1|2Q12I1 criteria. 
The solid line shows the approximate depth of the SWIRE 24/i.m 
imaging. 

lar mass. We use the observed 3.6^m flux to determine 
the rest-frame i^-band lumin osity, adopting an 1 1 Gyr 
single stellar population from iBruzual fc Charlotl (1200311 
to calculate the it'-correction. We then take the average 
Jf-b and mass-to4ight ratio of a red galaxy (M/L^ = 1) 
from iBell et al.l (j2003[ l. with a scatter of O.ldex. To de¬ 
termine the sSFR we furth er scale the masses by 1.65 to 
convert to a Salpeter IMF (jRaue fc Meved[20T^ . 

Figure [8] shows the sSFRs of the individual BCGs 
with redshift. This figure also contains the evolution of 
the sSFR of the m ain sequence of star forming galaxies 
(jElbaz et al.ll2011ll (solid line), as well as the division be¬ 
tween main sequence and starburst systems as defined 
by Elbaz et al. 

In Figure [9] we show the sSFR of the 24^m-detected 
BCGs as a function of stellar mass. The size of the points 
increases with the redshift of the BCG so that the largest 
points correspond to at z = 1.8 and the smallest points 
to z = 0.1. Overlaid we show the location of the main 
sequence of star forming galaxies in the field for the three 
redshift regimes that encompass our sample: z = 0,1, 2 
(jElbaz et al.ll2007t iDaddi et al.ll2007t l. again corrected to 
a Salpeter IMF. 

In both Figures Eland [9] general agreement is seen be¬ 
tween the sSFRs of the BCGs and the overall level seen in 
the field at a given redshift or mass. This indicates that 
although the star formation activity within the BCGs is 
quantitatively high, it is not in great excess of that seen 
for field galaxies, and indeed the BCGs may be classified 
as main sequence star forming systems. 

5. DISCUSSION 

5.1. A Change in the Observed Activity in BCGs with 

Redshift 

These data show a measurable change in the IR prop¬ 
erties of similarly selected BCGs from redshift z ^ 0.3 to 
z ^ 1.8. The IR colors indicate that the bulk of this evo¬ 
lution is due to an increase in the star-formation rates of 
BCGs to higher redshift, and the remainder from dusty 
AGN. This is qualitatively similar to the rapid increase 
of the global star formation rate density of the universe 
to z ^ 1, as well as the migration of to higher lumi- 



redshift 

Fig. 8. — The specific star formation rate of the 24/rm-detected 
BCGs shown as a function of r edshi ft. Overlaid (solid line) is the 
best fit relation of lElbaz et al.l (|2Qlin for main sequence star form¬ 
ing galaxies. The dotted line shows the rough division between 
starburst galaxies and main sequence galaxies, as classified by El¬ 
baz et al. No AGN correction has been made to the SFRs. 



Fig. 9. — The specific star formation rate of the 24/r.m-detected 
BCGs shown as a function of stellar mass. Overlaid are the re¬ 
lations fo r main sequence s tar forming galaxies at the r^ shifts 
z = 0,1 ilElbaz et al.|[200^ and z = 2 JPaddi et ahlf^OOTI^ that 
encompass our sample. The size of the data point, for each BCG, 
increases with increasing redshift, from z 1-1.8. 

nosities over the same redshift range. Figure E] suggests, 
however, that the BCGs may follow a steeper evolution 
than seen for the held. We cannot constrain the slope 
of this evolution, because we are only sensitive to the 
most luminous galaxies at high redshifts, however we do 
see clear increase in the fraction of BCGs detected above 
an IR-depth which is dehned relative to the evolution of 
the held IR luminosity function (see Section 4.2). Note 
that this result is not driven by the increase in AGN con¬ 
tamination: it remains when all of the candidate AGN 
are removed from the sample. This is evidence that the 
physics driving the increased activity in BCGs at higher 
redshift is either different than that driving the held evo¬ 
lution or is accelerated in the cores of galaxy clusters. 

5.2. The Stellar Mass Growth of BCGs 

Below z ^ I the SFRs of individual detected BGGs do 
not exceed 100 MQyr“^. If the detection rate of BGGs 







































































at this level (~ 10 %) is reflective of the duty cycle of the 
star-burst then they must be short-lived, or < 1 Gyr. 
At such rates, this star-formation will contribute a rela¬ 
tively small amount of stellar mass to the overall system. 
There have been conflicting measurements of the actual 
stellar mass growth rate of BCGs in this redshift range, 
but ge nerally an increase between ^ 0-2 x is seen since 
z = 1 (|Lidnian et al.l[20T^ iLiu et al.ll2012ll . This mass 
growth has been attributed by many to the accretion 
of established stellar populati ons through gas-poor ma¬ 
jor or minor galaxy mergers (lEdwards fc PattCT 120121 : 
iLidman et al.l l^lSll . The IR results seen here (below 
z < 1 ) do not add enough stellar mass to contradict this 
scenario. 

By z ~ 1, however, the star formation level has dras¬ 
tically increased to ^ 500-2000 M 0 /yr. Using similar 
arguments as above we reach similar timescale conclu¬ 
sions: if the duty-cycle is 20 % over the epoch studied, 
then the bursts must be limited to a few hundred million 
years. In this case however, the greatly increased star 
formation rates mean that even over this short timescale 
a BGG can easily double its stellar mass. At these levels 
in situ star formation is an important, perh aps dominant, 
contributor to the mass growth of BCGs (|Gollins et al.l 
120091: iStott et 311120111 : ILidman et al.ll201^ . 

This general reasoning can also be illustrated in Figure 
m which is adapted from ILidman et ^ (I2012D . Using 
NIR-determined mass measurements of BCGs Lidman et 
al. calculated the fraction of z = 0 stellar mass in BCGs 
in place at a given redshift, to z ~ 1. They find rela¬ 
tively good agreement with the build up of mass due t o 
dry-mergers in the SAMS of iDe Lucia fc Blaizod (1200711 . 
These two studies are shown in Figure [TOl but are plot¬ 
ted as the amount of mass still to be added to a BCG at 
a given redshift, as a fraction of the final mass. We com¬ 
pare to this the amount of stellar mass added by the star 
formation seen here, again in terms of the final stellar 
mass at z = 0 which is shown by the grey area. 

To determine the grey area we measure the average 
SFR for the entire BCG population in redshift bins, as¬ 
signing a SFR of 0 M 0 yr“^ to the undetected BCGs. 
The chosen bins are z =0.2-0. 6 , 0.6-1.0, 1.0-1.4, 1.4-1. 8 . 
We then integrate this star formation over the duration 
of the redshift bin to obtain the average amount of stel¬ 
lar mass added in that bin. This is then combined with 
existing average stellar mass of the BCGs within the red¬ 
shift bin to provide a total, final mass, treating each bin 
independently. This is complementary to our reasoning 
at the very beginning of this section where we took the 
high star-formation-rates of the IR-detected BCGs and 
used the detected fraction to constrain the limits of in¬ 
tegration (the duty cycle). Here we do not constrain the 
timescale of the star-formation -- allowing the star forma¬ 
tion to proceed at the same rate over the entire redshift 
bin but employ a lower average star formation rate 
which incorporates the IR-faint galaxies as well. The 
upper bound of the area includes all IR galaxies while 
the lower bound was computed with the candidate AGN 
removed. 

Figure 10 is meant only to illustrate the approximate 
importance of the dust enshrouded star formation seen 
here, relative to the amount of mass required by other 
observations. The stellar mass determination is crude 
and no attempt is made to correct for progenitor bias, 



Fig. 10.— This plot, which is adapted from ILidman et al.l II2012I ). 
shows the amount of stellar mass that is added to BCGs with 
time, normalized to their present-day {z = 0) stellar mass. Solid 
points denote the stellar mass measurements of BCGs from Lid- 
man et al., and th e solid line denotes the model predictions of 
IDe Lucia fc BlaizoH 11200711 . This is a modification of the Lidman 
et al. plot; here a value of 0.2 indicates that a BCG has a deficit of 
20% in its stellar mass, compared to today. The solid grey region 
shows the amount of stellar mass which, given the assumptions 
outlined in the text, would be added to BCGs by the dusty en¬ 
shrouded star formation seen in this work. The upper boundary is 
calculated by assuming all of the IR-bright systems detected here 
are star forming, while for the lower bound we have removed all 
possible AGN (Fig. |^ from the analysis. The intention of this plot 
is to compare the approximate mass growth through star formation 
with that required by the directly observed mass change (Lidman) 
and that expected from the models (De Lucia & Blaizot). 

as done in Lidman et al. Nevertheless, a simple picture 
emerges. Out to z ^ 1 in situ star formation adds less to 
the total stellar mass of the BCGs (as a population) than 
do the dry-merger predictions, and is less than required 
by the mass measurements of Lidman et al. Above z 1 
however, the situation is less clear. The amount of mass 
added by star formation rises steeply and begins to ap¬ 
proach the mass assembled through dry mergers in the 
models. At z ^ 1.5 both processes would more than dou¬ 
ble the existing mass of the BCG. More mass measure¬ 
ments of BCGs above z ~ 1 are required to determine if 
there is actual tension between these two models of mass 
assembly. Clearly both processes could operate and play 
significant roles in the formation of BCGs; however their 
combined effect is constrained by the observed buildup 
of stellar mass. 

5.3. Comparison to Recent Literature 

The results presented here stem from the first infrared 
analysis of a large sample of BCGs at 0.2 < z < 1.8 
They are, however, in solid agreement with complemen- 
tary studies of ot h er ma ssive galaxies and single clusters. 
iMarchesini et al.l (120141) undertook a study of the pro¬ 
genitors of local Ultra Massive Galaxies (UMGs) within 
the UltraVista Survey using abundance matching tech¬ 
niques. While not explicitly identified with clusters, the 
stellar-mass of 10 ^^ ® Mq is similar to the mass of our 
BCGs 3x10^° - 7x10^^, and it is likely the Marchesini et 
al. population overlaps significantly with the population 
studied here. They find remarkably consistent results: 
the progenitors of UMGs remain quiescent until z ^1, 
beyond which the exhibit extreme levels of dusty star 
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formation of several hundreds of solar masses per year. 

Extreme star formation activity has also been mea¬ 
sured tow ard or in th e core s of several high-redshift 
clusters. iSantos et al.l ()2ni5[ l reports a total SFR of 
1875M0yr“^ within the central 250kpc of a z =1.58 
cluster, distributed over three massive galaxies. Indeed, 
many groups are now finding the star formation activ¬ 
ity within the central r egions of z >1.5 c l usters, though 
not a.lwavs the BCG jiTran et al.1 l201fll iPooesso et al.l 

120121 : iFassbender et al.l[^llLl20l4lBrodwin et al.ll2013ll . 

is co nsistent with the level o f field activity. In some cases 
(e.g. IFassbender et al.ll20l4l this activity is due to an en¬ 
hancement in galaxy-galaxy mergers within the cluster 
cores. 

5.4. Star Formation Driver: Dissipative Mergers? 

These results of the previous sections indicate a rapid 
increase in star formation in the central galaxies of clus¬ 
ters beyond z ^ 1 and this in turn suggests a increase 
in the efficiency of gas deposition onto these systems. In 
the field, major galaxy mergers appear to be the domi¬ 
nant method of delivering new (and large) reservoirs of 
gas to galaxies, but this process may be problematic in 
very high density regions. Galaxies within the cluster en¬ 
vironment are subject to ISM removal processes such as 
ram pressure stripping or strangu lation and may be gas 
deficient le.g. lGavazzi et al.ll2006f) by the time t hey re ach 
the centre. As pointed out in lMcDonald et ^ (I2012fl for 
the Phoenix Gluster BCG, the prodigious star formation 
rates measured here require the accretion galaxy with a 
gas mass far in excess (lO-lOOx) that of a typical (z = 0) 
late-type cluster galaxy. Certainly in-falling field galax¬ 
ies at high redshift will begin their accret ion with more 
gas than at low redshift (jGeach et al.ll20lTh and therefore 
may be able to retain a larger overall mass of gas during 
their infall. Indeed, simulations of galaxy-ICM interac¬ 
tion indicate galaxies may retain substant ial fractions of 
their gas during infall to the cluster center (|Jachvm et al.l 
1200711 , but the gas content of cluster galaxies beyond z ^ 
1 remains unconstrained. 

Still, the one high-redshift and 24yLtm-bright for which 
we have conducted high spatial resolution i maging with 
the H ST shows clear indications of a merger (iWebb et al.l 
1201511 . This system, SpARCS1049-|-56, consists of a 
long ( 66 kpc) tidal tail interspersed with clumps along 
its length. It appears to originate from within the stellar 
halo of the BCG and exhibits copious amounts of star 
formation (^1100 MQyr~^). This gas-rich merger would 
not have been apparent without HST NIR imaging and 
thus ruling out this scenario for the rest of the sample 
is premature. If SpARCS1049-|-56 is indeed representa¬ 
tive of the z >1 24/rm-bright SpARCS sample, then wet 
gas-rich mergers can occur at the centres of clusters and 
indeed, are an important process in building the stellar 
mass of BCG at early times. 

On the other hand, studies of cool-core clusters 
(|Rawle et al.ll2012t iLiu et al1l2012tl reveal a strong corre¬ 
lation between the star formation rates of BCGs (limited 
to below z ~ 0.5) and their cooling times (as measured 
in the X-ray). The sample presented here is optical- rich¬ 
ness selected (nominally stellar-mass-selected) and there 
is no information on the state of the gas in these sys- 
tems. X-ray studi es of the Red-Sequence Cluster Survey 
((Hicks et al.l 120131) . which are limited to z < 1, indicate 


that perhaps 10 % of optically selected clusters harbour 
cool-cores. This fraction is enough to account for all 
IR-detected galaxies in our sample to z ^ 1, but does 
require the number of cool-core clusters to increase sub¬ 
stantially at higher redshifts. A detailed comparison of 
SpARCS optical and X-ray selected clusters has not yet 
been made at these redshifts. One SpARCS field (XMM- 
LSS) has X-ray coverage and we have compared the clus- 
te rs found by the SpA RCS survey with those identified 
in lAdami et al.l ((201 If) and IWillis et al.l (I2013f) . At the 
significance and richness limit of SpARCs the overlap is 
primarily restricted to the most X-ray luminous systems 
at z < 0.6, where the IR detection rate is low. Only one 
common cluster has a BCG detected at 24^m and it is 
not unusual in its X-ray luminosity. X-ray observations 
of the SpARCS sample would provide a simple yet clear 
confirmation of this idea - the 20 % of the sample which is 
IR-luminous should be cool-core clusters when observed 
in the X-ray. 

6. SUMMARY, CONCLUSIONS AND FUTURE WORK 

We investigate the infrared properties of a large sam¬ 
ple (535) Brightest Cluster Galaxies in massive opti¬ 
cally selected galaxy clusters over the redshift range 
0.2 < z < 1.8, using the SWIRE public 24/xm imaging 
and catalogs. We find several important results: 

• We detect 125/535, or 23% of the BCGs at 24/rm. 
The fraction of BCGs detected at 24/xm above 
100/rJy and above a fixed IR luminosity increases 
beyond redshift z ^ 1. Below z < 1 no BCGs show 
Lir > 10 ^^ Lq, whereas above z > 1 this fraction 
rises to 20 %. 

• An investigation of the Spitzer-IRAC IR colors in¬ 
dicates that the bulk of the IR emission of BCGs is 
not AGN dominated, but rather is due to dust en¬ 
shrouded star formation. Only 7/125 of the 24/im- 
detected BCG sample inhabit the AGN region (de¬ 
fined by Donley et al.(2012)) of the IRAC color- 
color diagram. 

• Assuming no contribution from AGN to the in¬ 
frared luminosity, we calculate the star-formation- 
rates for the BCGs. These range from ~1 M 0 yr“^ 
at low redshift to ^ lOOOM 0 yr —1 at high redshift. 
At all redshifts, however, the specific SFRs of the 
BCGs are consistent with those of main-sequence 
star-forming galaxies of similar mass, and thus the 
bulk of the BCG population is not starbursting. 

• We argue the the star formation episodes are short¬ 
lived and below z ~ 1 do not contribute more than 
10% to the final total stellar mass of the BCGs. In 
this case BCG growth is likely dominated by dry 
mergers. Above z ~ 1 however, the star formation 
seen here may double the stellar mass of the BCGs 
and is thus an important, and perhaps dominate, 
process for stellar mass growth. 

• The physics driving the increase in star formation 
above is as yet unconstrained, although the one ob¬ 
ject for which we have deep and high spatial resolu¬ 
tion imaging (SpARCS1049-|-56) shows clear signs 
of wet merger activity. 
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TABLE 2 

Infrared Detected BCGs 


RA 

Dec 

redshift® 

S24/4m (l4lJy) 

Lir (loll L0)b 

SFR (Meyr-l) 

Stellar Mass (10^^ M©) 

10:47:22.8 

57:00:51 

0.238 

1500.4±15.5 

0.71±0.36 

12±6 

0.9 

10:33:34.0 

58:14:35 

0.430 

17900.7±28.6 

16.1±1.2 

281±218 

1.0 

10:49:22.3 

56:56:43 

0.244 

1429.9±16.7 

1.9±1.3 

33±24 

2.4 

10:48:50.3 

56:09:08 

0.468 

187.0±17.3 

0.34±0.24 

6±4 

1.2 

10:42:24.6 

57:59:49 

0.387 

627.3±16.5 

1.1±0.8 

19±13 

0.6 

10:39:30.1 

57:14:35 

0.425 

386.9±18.6 

0.94±0.58 

17±10 

0.6 

10:50:56.6 

57:07:20 

0.475 

221.4±12.2 

0.66±0.40 

12±7 

0.6 

10:37:38.9 

58:45:29 

0.688 

329.0±16.8 

2.2±0.7 

39±12 

1.8 

10:51:57.3 

56:59:05 

0.700 

638.7±12.8 

4.4±1.4 

77±24 

1.0 

10:44:48.9 

59:12:16 

0.737 

239.0±12.7 

1.9±0.5 

32±8 

1.0 

10:39:46.1 

58:54:17 

0.750 

1380.2±17.1 

11.0±0.4 

192±72 

2.8 

10:42:28.6 

57:51:24 

0.778 

464.0±17.6 

4.1±1.1 

71±19 

1.0 

10:37:25.8 

59:02:12 

0.813 

248.5±13.5 

2.3±0.6 

39±10 

1.4 

10:42:36.0 

58:18:16 

0.825 

163.3±14.8 

1.5±0.4 

26±7 

1.0 

10:40:05.0 

59:16:26 

0.825 

266.7±16.9 

2.5±0.7 

43±12 

1.4 

10:35:36.5 

57:50:35 

0.875 

175.4±15.8 

1.7±0.6 

29±10 

1.1 

10:58:24.5 

57:27:26 

0.888 

234.1±15.5 

2.3±0.8 

41±14 

1.3 

10:46:43.4 

58:21:20 

0.900 

411.5±13.9 

4.5±1.5 

78±26 

2.4 

10:48:17.1 

55:52:18 

0.900 

785.9±16.5 

9.1±3.5 

159±62 

1.0 

10:44:05.0 

58:59:04 

0.925 

514.7±11.9 

6.2±2.2 

109±39 

1.4 

10:51:55.3 

57:01:09 

0.937 

269.5±11.2 

3.2±1.2 

56±20 

0.6 

10:55:56.3 

57:10:37 

0.950 

560.9±17.5 

7.4±2.9 

130±50 

0.8 

10:34:05.4 

58:11:31 

0.963 

523.4±18.1 

7.2±2.9 

126±51 

1.9 

10:52:22.7 

56:55:44 

0.963 

525.8±11.6 

7.3±2.8 

127±49 

1.1 

10:46:25.0 

57:04:29 

1.000 

189.6±15.2 

2.7±1.3 

47±22 

1.7 

10:47:05.7 

58:05:42 

1.013 

2551.0±19.7 

54.1±20.3 

947±356 

2.6 

10:52:06.0 

56:56:30 

1.038 

485.5±11.8 

8.8±3.8 

153±66 

1.7 

10:47:49.7 

57:43:37 

1.029 

388.9±14.1 

6.9±2.8 

120±50 

1.0 

10:47:38.8 

57:52:29 

1.075 

1281.6±18.0 

32.5±12.2 

570±213 

2.0 

10:59:55.7 

57:45:56 

1.075 

431.0±17.9 

8.9±3.9 

155±69 

1.2 

10:52:30.0 

58:19:43 

1.125 

357.4±15.5 

8.6±3.7 

150±65 

1.6 

10:46:26.4 

57:32:30 

1.125 

926.6±16.0 

27.0±10.1 

472±177 

1.5 

10:59:14.9 

57:32:40 

1.200 

226.4±14.4 

6.6±3.1 

116±54 

1.4 

10:52:08.4 

57:02:38 

1.250 

377.2±11.1 

15.6±0.7 

272±123 

2.2 

10:50:07.2 

57:16:51 

1.600 

662.7±12.3 

82.2±24.3 

1438±426 

3.6 

10:44:59.5 

57:52:07 

1.650 

281.1±17.0 

34.2±9.0 

599±152 

5.3 

10:49:22.6 

56:40:33 

1.70 

606.1±18.0 

104±29 

1819±511 

4.5 

02:19:41.8 

-04:00:33.3 

0.141 

358.1±18.4 

0.12±0.06 

2±1 

0.5 

02:25:24.9 

-03:47:35.4 

0.250 

263.8±16.3 

0.18±0.09 

3±2 

0.3 

02:24:26.9 

-05:36:32.2 

0.262 

1814.7±17.8 

1.1±0.5 

19±9 

0.7 

02:24:33.1 

-04:53:56.1 

0.337 

262.34±19.7 

0.34±0.29 

6±5 

0.5 

02:19:10.5 

-03:43:34.4 

0.763 

2684.4±18.8 

18.8±7.6 

329±132 

1.6 

02:20:13.3 

-06:00:54.9 

0.688 

339.5±20.9 

2.3±0.8 

40±13 

1.1 

02:20:27.8 

-05:47:24.9 

0.762 

236.0±16.4 

1.9±0.5 

34±9 

1.2 

02:22:13.0 

-04:21:58.4 

0.788 

318.1±20.2 

2.8±0.8 

48±13 

0.9 

02:24:24.0 

-02:58:02.7 

0.788 

216.3±19.6 

1.9±0.5 

33±9 

0.9 

02:19:29.0 

-04:07:00.2 

0.813 

1086.4±20.4 

10.3±3.8 

180±66 

1.6 

02:18:34.4 

-05:00:43.6 

0.651 

368.8±19.0 

3.5±1.0 

61±17 

1.7 

02:15:43.9 

-04:24:53.5 

0.850 

694.6±18.2 

7.±2.0 

122±37 

2.7 

02:24:29.2 

-04:10:13.1 

0.900 

1082.9±19.8 

12.4±5.4 

235±95 

1.3 

02:25:06.9 

-04:47:18.4 

0.925 

276.0±19.4 

3.2±1.2 

55±21 

1.8 

02:02:08.5 

-03:41:26.8 

0.925 

345.9±19.3 

4.0±1.5 

70±27 

1.0 

02:14:38.2 

-03:37:38.3 

0.931 

2241.6±22.2 

35.1±12.8 

614±224 

2.7 

02:23:05.8 

-04:13:35.5 

1.048 

240.1±18.2 

3.1±1.3 

53±23 

0.9 

02:20:54.8 

-03:32:57.5 

0.994 

1137.5±21.7 

20.9±8.8 

356±153 

3.3 

02:22:36.5 

-03:50:30.3 

0.822 

1039.2±19.7 

19.3±8.4 

338±146 

2.0 

02:16:38.2 

-03:28:44.6 

1.00 

698.5±21.1 

11.3±5.5 

198±95 

2.0 

02:22:54.3 

-04:14:11.9 

1.00 

362.3±21.3 

5.5±2.4 

96±42 

1.0 

02:16:03.2 

-03:33:57.9 

1.05 

352.5±16.4 

6.5±2.7 

113±47 

1.8 

02:18:05.2 

-05:00:10.5 

1.095 

247.1±15.9 

4.8±2.1 

83±36 

2.2 

02:16:52.3 

-03:37:56.7 

1.350 

508.8±20.0 

32.0±11.7 

559±204 

2.3 

02:26:15.6 

-04:56:28.0 

1.375 

182.9±19.7 

10.3±4.9 

180±87 

2.3 

02:21:43.8 

-03:21:57.0 

1.375 

1023.9±21.0 

72.2±27.1 

1263±474 

3.8 

02:27:30.5 

-04:32:03.4 

1.450 

435.9±18.2 

35.8±11.1 

628±195 

1.8 

02:18:45.2 

-05:42:56.8 

1.45 

313.3±19.6 

24.8±8.3 

434±145 

2.0 

16:05:44.6 

54:57:16.2 

0.188 

378.6±16.8 

0.01±0.06 

2±1 

0.7 

16:06:27.9 

54:56:14.7 

0.188 

231.0±12.9 

0.08±0.04 

1±0.7 

0.3 

16:13:32.9 

56:17:47.5 

0.212 

683.5±17.4 

0.27±0.13 

5±2 

0.8 

16:18:38.3 

55:17:13.7 

0.250 

107.5±13.2 

0.08±0.04 

1±1 

0.9 

16:14:51.5 

54: 02: 22.0 

0.275 

1042.6±16.7 

0.72±0.36 

13±6 

0.5 

16:19:14.8 

55:13:40.8 

0.313 

224.6±16.4 

0.21±0.21 

5±4 

0.7 

16:10:11.2 

53:46:24.5 

0.362 

290.0±15.9 

0.44±0.37 

8±7 

0.9 

16:13:17.5 

56:01:19.3 

0.387 

2571.9±18.8 

3.9±2.8 

68±49 

1.3 

16:00:37.2 

55:27:26.9 

0.465 

163.5±16.2 

0.33±0.23 

6±4 

1.4 

16:08:57.2 

56:00:18.6 

0.688 

524.1±17.0 

3.5±1.1 

61±20 

7.2 



TABLE 2 


Continued 


RA 

Dec 

redshift® 

S24/im (mJy) 

Lir (loll L0)b 

SFR (Mgyr-l) 

Stellar Mass (10^^ M©) 

16:02:32.3 

54:56:59.6 

0.700 

150.0±14.5 

1.1±0.3 

18±6 

0.8 

16:18:40.1 

54:50:48.0 

0.700 

759.9±15.1 

5.2±1.7 

91±29 

1.2 

16:06:28.8 

55:33:07.7 

0.725 

165.3±16.1 

1.2±0.3 

22±6 

0.9 

16:18:26.4 

54:58:28.7 

0.737 

278.2±13.0 

2.1±0.6 

38±10 

1.3 

16:11:01.3 

54:17:05.3 

0.750 

695.7±15.7 

5.5±1.6 

97±28 

1.8 

16:03:04.5 

54:57:24.3 

0.775 

607.1±15.9 

5.2±1.4 

90±25 

1.2 

16:09:29.4 

54:29:40.6 

0.788 

1764.0±16.7 

17.3±6.0 

304±106 

0.8 

16:17:29.9 

55:58:18.2 

0.788 

207.1±16.2 

1.8±0.5 

31±8 

1.0 

16:13:05.7 

55:59:50.0 

0.825 

169.9±13.8 

1.5±0.4 

27±7 

1.4 

16:07:16.8 

55:32:59.8 

0.825 

600.2±13.8 

5.7±1.6 

100±28 

1.1 

16:08:25.6 

54:45:08.9 

0.906 

610.6 ±13.4 

5.8±1.6 

102±29 

0.8 

16:19:19.7 

55:36:45.7 

0.862 

299.8±13.9 

2.9±0.9 

51±15 

0.6 

16:02:50.7 

54:54:52.2 

0.862 

337.9±14.7 

3.3±1.0 

58±17 

0.9 

16:12:10.4 

56:07:53.6 

0.888 

483.1±13.9 

5.2±1.7 

90±29 

0.8 

16:11:12.7 

55:08:23.6 

0.907 

3065.4±19.9 

43.7±15.8 

766±275 

1.8 

16:08:07.8 

54:19:43.1 

0.900 

765.5±17.5 

8.8±3.3 

154±57 

1.4 

16:18:18.0 

54:35:38.9 

0.944 

671.7±17.0 

8.8±3.7 

155±66 

1.8 

16:22:42.1 

54:50:55.0 

1.000 

390.7±19.7 

6.0±2.6 

105±45 

1.1 

16:02:43.5 

54:49:09.7 

1.013 

194.3±17.4 

2.9±1.4 

51±24 

0.8 

16:13:01.7 

54:46:10.0 

1.092 

1220.8±14.5 

24.2±9.5 

423±167 

1.8 

16:08:27.1 

54:36:47.3 

1.125 

345.3±14.0 

8.3±3.3 

145±58 

1.2 

16:07:25.7 

54:40:41.2 

1.125 

250.6±15.0 

5.8±2.5 

101±43 

0.7 

16:11:15.8 

54:15:11.4 

1.125 

365.7±15.6 

8.8±3.9 

153±69 

1.3 

16:05:55.2 

55:31:47.2 

1.350 

326.0±14.2 

19.5±8.0 

342±138 

2.6 

16:06:52.7 

55:39:36.6 

1.650 

241.9±16.5 

28.3±7.3 

496±128 

4.6 

16:40:33.2 

41:23:11.7 

0.662 

220.6±16.2 

1.4±0.5 

25±8 

1.6 

16:34:35.3 

41:36:14.6 

0.781 

584.5±15.1 

4.2±1.3 

73±22 

1.3 

16:38:22.9 

41:53:19.1 

0.862 

772.1±16.4 

8.0±2.5 

141±44 

1.3 

16: 37:35.1 

40:17:18.1 

0.912 

217.1±16.3 

2.3±0.9 

40±16 

2.9 

16:41:52.4 

41:34:01.5 

0.925 

297.8±17.8 

3.4±1.3 

60±23 

0.7 

16:37:49.5 

42:04:55.3 

0.931 

2767.4±18.5 

43.7±16.1 

765±281 

3.7 

16:37:01.1 

41:31:02.4 

0.944 

747.8±17.1 

10.0±4.5 

175±78 

1.3 

16:37:46.5 

40:37:32.3 

0.944 

507.1±17.4 

6.5±2.5 

114±44 

2.3 

16:31:06.1 

40:50:31.3 

0.956 

340.8±17.8 

4.3±1.8 

76±31 

1.1 

16:38:07.1 

41:57:34.2 

1.075 

362.9±16.4 

7.3±3.0 

128±53 

0.8 

16:40:22.3 

42:03:14.4 

1.125 

295.8±13.3 

7.0±2.9 

122±50 

1.2 

16:37:27.0 

41:57:00.6 

1.125 

334.7±13.3 

8.0±3.2 

139±57 

1.3 

16:41:40.3 

41:13:18.4 

1.125 

283.4±15.0 

6.4±2.8 

116±49 

1.3 

16:39:12.2 

40:57:25.4 

1.275 

269.8±16.1 

10.9±5.7 

190±100 

1.7 

16:37:59.2 

42:03:56.5 

1.575 

175.0±12.9 

17.0±5.0 

293±87 

1.7 

16:37:37.9 

41:17:28.6 

1.775 

339.4±17.9 

52.9±14.3 

925±249 

2.2 


Bold denotes a spectroscopic redshift, non-bold corresponds to the red sequence estimated redshift. 
LIRs estimated using the 24/im flux and ICharv ElbazI l|2Q01il methodology. 











